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Abstract 

This paper addresses two aspects of duct 
propagation and radiation which can contribute to 
more efficient fan noise predictions. First, we assess 
the effectiveness of Rayleigh’s formula as a ducted But 
noise prediction tool. This classical result which 
predicts the sound produced by a piston in a flanged 
duct is expanded to include the uniform axial inflow 
case. Radiation patterns using Rayleigh’s formula with 
single radial mode input are compared to those 
obtained from the more precise ducted flan noise 
prediction code TBIEM3D. Agreement between the 
two methods is excellent in the peak noise regions both 
forward and aft. Next, we use TBIEM3D to calculate 
generalized radiation impedances and power 
transmission coefficients. These quantities are 
computed for a wide range of operating parameters. 
Results were obtained for higher Mach numbers, 
frequencies, and circumferential mode orders than 
have been previously published. Viewed as functions 
of frequency, calculated trends in lower order inlet 
impedances and power transmission coefficients are in 
agreement with known results. The relationships are 
more oscillatory for higher order modes and higher 
Mach numbers. 
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m = 0£1£2 „ . . Circumferential mode number 

M Inflow Mach number 

n Radial mode number 

r D Duct radius 

t Time 

p' acoustic pressure in ground fixed frame 

P m m-th modal coefficient of acoustic pressure 

in duct fixed frame 

(r.yr.Z) Cylindrical coordinates in duct fixed frame 


U " m-th modal coefficient of axial component of 

acoustic velocity in duct fixed frame 

v 0 normal component of acoustic velocity in 
ground fixed frame 

P - Jl-M 1 2 Compressibility (stretching) 

parameter 
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v m n-th zero of J' m , n= 0,1,.. . for m = 0 and 
n = 1 . 2 ,... for m * 0 

IT" Power transmission coefficient for m-th 

circumferential mode and n-th radial mode 

p 0 Ambient density 

Generalized radiation impedance for m-th 
circumferential mode due to interaction of n- 
th radial pressure mode with l-th radial 
velocity mode 

Q Temporal frequency (radians/second) 

Introduction 

Today's civil jet transport aircraft are about 20 dB 
quieter than were those introduced some 40 years ago. 
Aircraft became quieter as turbojets were replaced by 
turbofans, sound-absorbing liners were added to the 
inlet and exhaust ducts, inlet guide vanes were 
eliminated, and bypass ratios grew larger. Still, airport 
communities around the world would like to see even 
less noisy aircraft. 

As bypass ratios have increased, average exhaust 
velocities have decreased, resulting in less jet noise. At 
the same time, fan noise has increased because more of 
the thrust comes from the fan. 

This paper addresses two aspects of duct 
propagation and radiation which can contribute to 
more efficient fan noise predictions, facilitating 
parametric screening and trend studies. In both cases, 
we use the ducted fan noise prediction computer 
program TBIEM3D 1 . The code, based on a boundary 
integral equation method (BIEM) 2 , calculates the 
acoustic pressure anywhere in the sound field due to 
the scattering of incident sound by a finite length, 
cylindrical duct in a uniform flow field. The duct walls 
can be hard or lined and incident sound is generated by 
simple point sources. 

First, we assess the effectiveness of Rayleigh’s 
formula as a ducted fan noise prediction tool. The 
formula gives the acoustic pressure field due to a piston 
in a rigid wall and requires the acoustic pressure or the 
normal component of velocity on the piston as input. 
For the no inflow case and when only one propagating 
mode is present the formula is simple to implement 
and has been applied to inlet radiation problems by 
many researchers. 

In this paper, we extend the formula to include the 
uniform axial inflow case. Noise radiation results 


using Rayleigh's formula widi single propagating 
mode input are compared with TO I EM 3D calculations. 

Next, we use TBIEM3D to calculate generalized 
radiation impedances and power transmission 
coefficients. When a propagating spinning mode in a 
duct reaches the inlet, part of the energy is radiated and 
part reflected due to the abrupt change in geometry. 
The process is complicated by edge diffraction. 
TBIEM3D accounts for these phenomena. Previous 
researchers have relied on less exact noise prediction 
methods 3 ' 8 for impedance calculations. Comparisons 
with these works are discussed. 

Generalized radiation impedances and power 
transmission coefficients provide important noise 
radiation characteristics to the engine designer. We 
demonstrate the usefulness and versatility of 
TBIEM3D as a design tool by calculating these 
parameters for a wide range of operating parameters 
and circumferential mode numbers. 

Part I: Rayleigh’s Formula 


Theory 


In the analysis that follows all variables have been 
nondimensionalized: length by r D , mass by p g rj , and 

time by r D c~’ . 

Rayleigh’s formula gives the acoustic pressure 
field due to a vibrating piston in a rigid wall 9 . Tyler 
and Soffin, in their classical work on duct propagation 
and radiation, proposed the formula for calculation of 
the acoustic radiation field due to propagating modes 
from a semi-infinite, flanged duct with no flow 10 . 

The application of Rayleigh’s formula to inlet (or 
exhaust) radiation requires knowledge of either the 
normal component of acoustic velocity or the acoustic 
pressure on the inlet (or exhaust) surface. It is 
assumed that these quantities are zero everywhere else 
on the plane containing the inlet (or exhaust). 
Furthermore, the acoustic pressure is assumed to be 
zero in the region of space behind the inlet disc if inlet 
radiation is considered or in front of the exhaust disc if 
exhaust radiation is desired. 

We briefly derive a version of Rayleigh’s formula 
that is valid with or without flow. Consider the 
radiation of sound from the inlet of an infinitesimally 
thin, hard walled, finite length, cylindrical duct in the 
presence of a uniform axial inflow. We adopt the usual 
spinning mode representation and assume that the m-th 
circumferential mode of acoustic pressure has the form 


p m (r,tff,Z,t)= P m (r,Z)e 


i kt-my 




( 1 ) 
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It is shown in reference 2 that the modal 
coefficient P m is a solution of the two-dimensional 
Helmholtz equation 


i d ( d 

r 


r dr 


where 


d 2 m 2 

, + — r + 

dr J d 2 % r 




\p m (r,%) = 0, (2) 


(3) 


The free-space Green’s function for the operator 
(2) is given by 


1 r e ~ iKR 

G(r,r',Z-Z') = — jcosm\i/—d\ff , 


(4) 


where 


R = -J7 2 + r' 2 -2rr'c(wyr + (< • (5) 

Assume that the inlet plane is located at the axial 
position and that the modal functions P m 

d P m 

and — — are known on the inlet disc and are zero 

everywhere else on the inlet plane. Then using the 
Helmholtz integral representation, we can write the 
acoustic pressure anywhere in the sound field as 11 


i 

P m (r,Z) = -jr'^(r\$Mr,r'4-S ln )d r '± 


1 

I 


( 6 ) 


r’P n {r'.^)^{r,r'^-^)dr' 


The positive sign in (6) applies to | and the 
negative sign to 

Evaluation of (6) is simplified considerably if 
predictions for £ < are not desired. For this 
situation, we set 

P m {r^) = 0 4 < l n (7) 

in (6) yielding the result 

F M (r,|) = 

i 

2 jr'P m (r',4 ln )^( r .r',4-^)dr' $>&.' (8) 

o 

Thus, only the acoustic pressure at the inlet is required 
to compute the forward radiation field. 


Denote the axial coordinate of the duct trailing 
edge by and apply the same procedure to the 
exhaust plane, then we obtain the composite function 

P m (r4) = 

0 

1 

• 2jr'P m (r\4 bl )^(r,r'4-4 in )dr' £>£„ (9) 

0 

1 

-lj r'P”(r\U^(ry Z < 

0 

for the modal coefficient of acoustic pressure. 

Equation (9) extends the classical Rayleigh 
formula. It includes radiation from the inlet and 
exhaust planes and is valid for uniform axial inflow. 

In practice, Rayleigh’s formula is usually applied 
to a single radial mode obtained from infinite duct 
theory. We assess the accuracy of this noise prediction 
technique by comparing radiation patterns computed 
using (9) with single radial mode input to those 
generated by TBIEM3D. 

Resu lt s 

In all calculations the duct radius and length are 
one meter. Incident sound is generated by a circular 
array of twenty spinning axial dipole point sources. 
The source disc is perpendicular to the duct axis and 
located in the middle of the duct. This configuration 
simulates the thrust component of loading noise 
generated by a fan with twenty blades. The rotational 
speed of the source disc is chosen so that only one 
radial mode is cut-on. 

In this study, we compare the radiated field from 
TBIEM3D with that obtained from (9) with single 
radial mode input. Two cases are examined, one with 
no inflow and the other with axial inflow Mach 
number M = 0.4. The rotational speed and source 
strengths of the dipoles are kept constant for the two 
cases. In order to obtain a proper comparison between 
the two methods, the amplitude of the single radial 
mode must be adjusted to match the TBIEM3D results. 
This is accomplished by a procedure described in Part 
II of this paper. 

In figure 1 the no inflow case is considered. We 
plot contours of sound pressure level for the m = 20 
circumferential coefficient of acoustic pressure in the r- 
Z plane. Figure la shows the field calculated using 
Rayleigh’s formula. TBIEM3D results are displayed in 
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figure lb. The radiated patterns for the two cases are 
very similar in the peak noise region in both direction 
and amplitude. 

Figure 2 contains equivalent results for the M = 
0.4 case. Again the agreement between the methods is 
excellent in die peak noise region both forward and aft. 

Since edge effects are ignored, details of the 
complicated diffraction pattern observed in the 
TBIEM3D results are not captured by the Rayleigh 
formula predictions. 

Discussion 

Rayleigh’s formula is simple to evaluate and 
involves little computational effort to produce the 
radiated field. For the single propagating mode case, 
the formula provides an excellent noise prediction tool 
in the peak noise region provided both the reflected 
and transmitted amplitudes of the mode are known. 
These quantities are not always available, but for many 
situations reflection can be ignored and the field 
accurately calculated with just the transmitted portion 
of the radial mode. 

Since the input function for Rayleigh’s formula 
can have an arbitrary radial distribution, (9) can also 
be used if more than one mode is propagating. One 
must then know all modal amplitudes. 

TBIEM3D does not explicitly resort to the radial 
decomposition of the pressure field inside the duct. In 
fact, near the duct edges the actual pressure field 
cannot be written as a sum of radial modes in the 
conventional sense. This is because of the complicated 
diffraction pattern caused by the sharp duct edges. 
These effects are accounted for implicitly by TBIEM3D 
and ignored by the Rayleigh formula representation. 
Thus, if fine scale details of the radiated field are 
required, then TBIEM3D would be the prediction 
method of choice. 

Part II: Generalized Radiation Impedance 


Theory 

The calculation of generalized radiation 
impedances requires knowledge of the acoustic 
pressure and acoustic velocity at the inlet. In this 
work, these quantities have been calculated by 
TBIEM3D. 

The concept of generalized radiation impedance is 
based on classical infinite duct, spinning mode 
representations of the acoustic pressure and velocity 
fields inside the duct. The m-th circumferential mode 
of acoustic pressure inside the duct has the form 


>(*< 

p m {r,\ l /,Z,t)= P m {r.Z)e [ * \ (10a) 


where 


P m {r,Z) = 

+ A Le^ z )J m ( Vmi ry 


(10b) 


and A ^ are as yet undetermined constants. Note that 
the summation in (10b) begins at n = 1 if m * 0 . 

The axial component of acoustic velocity is related 
to the acoustic pressure via the acoustic momentum 
equation 


ikuS - M 


du£ 

dZ dz 


= 0 


( 11 ) 


Using (10-11) u ” can be written as 


Uz(r,\f/,Z,t) = U”(r,Z)e 0 1 (12a) 

with 

Uz(r,Z) = 

X {M~ m Kne~‘ KmZ + /C K«f iK ~ z )j m {v m r) (12b) 

n=0 


and 




PMk^+k ‘ 


(13) 


The expansion (10) is not valid near the open ends 
of a finite length duct. This is because of the singular 
behavior of the acoustic pressure at the leading and 
trailing edges 2 . However, the acoustic pressure can be 
adequately represented by (10) for regions of space 
inside the duct that are sufficiently far from the ends,. 

TBIEM3D does not explicitly calculate the 
expansion (10). The axial coefficients A^ can be 

determined by calculating the acoustic pressure on 
several discs inside the duct. Hankel transforms are 
calculated for each disc to determine the strengths of 
the radial modes. The axial coefficients are then 
obtained by applying a least squares process to the 
Hankel transform results. Finally the axial component 
of acoustic velocity inside the duct is determined from 
(12-13). 


Inlet Impedance 

For the m-th circumferential mode, we define inlet 
impedance, £ m , by the ratio of acoustic pressure to 
axial velocity evaluated at the inlet. It is evident from 
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(10) and (12) that is a function of the radial 
variable only. 

Assume that the inlet disc is located at Z =0 and 
define the quantities P” and U” by 

04) 

and 

. (15) 

These known quantities are related through the inlet 
impedance, i.e.. 


X [/T - U"C m (r)\j m {v m r) = 0 . (16) 


The generalized impedances, are obtained by 
applying the orthogonality relation 12 


i i 

}rJ m (v mi r K(v m » r )dr = - 
0 * 


( 


1 - 


m 


,(17) 


where 8 is the Kronecker delta, to (16). This yields 
the result 


write the acoustic pressure and axial velocity for a 
single radial mode as 


P m (r,Z ) = +R:e"*- Z )j m (v mn r ) , (22) 


and 




U”(r,Z) = 

+/C*;* +1,C - Z 




(23) 


where the pressure reflection coefficients, 


given by 



R* are 


(24) 


and account for reflection from the duct inlet. 

The inlet impedance for this situation is constant 
and given by 


r 


l + R” 


/*«+j uLK 


(25) 


Using (22-25), we write the axial intensity at the inlet 
in terms of the inlet impedance. Thus, 


pm 
* n 



(18) 


where the generalized radiation impedances are given 
by the formula 



/ 


\rS m {r)j m ( Vml r)j m {v^r)dr 


I z\ inlet =4l^ / «( V - r f 1^- >< 

1 t . (26) 

Re{{C m - Mjl-MC)} 

We are also interested in the axial intensity at the 
inlet due to the incident wave only which we denote by 
(/”) | w „ . This is obtained from (26) by setting 

R" = 0. Therefore, 


Acoustic Power Transmission Coefficient 

The acoustic power radiated from the inlet due to 
the m-th circumferential mode is obtained by 
integrating the axial component of the acoustic 
intensity vector over the inlet disc. 

At any location inside the duct, the axial 
component of intensity can be written as 

l” =-j/?e{(p w - MU?)(U2 - AfP m )*} . (20) 

In (20), the asterisk superscript denotes the complex 
conjugate. Evaluating the axial intensity at the inlet 
yields 

=- L 2 \Vz( MC)'} • ( 21 ) 

The power radiation analysis that follows will be 
considered for individual radial modes. Using the 
axial coefficients obtained from TBIEM3D above, we 


teLU = 4 | x 

2 ^ . (27) 

Define the acoustic power transmission coefficient 
for the (m,n) mode, n* , as the ratio of power 
transmitted through the inlet disc to incident power 
through the disc. For the case considered here, this 
quantity is simply the ratio of axial intensities at the 
inlet which yields the result 

n;=|p»+tfjcfx 

Re{{Z m -M\l-Mt; n )‘} ' . (28) 

Re^l-Mn^^-M)') 
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1‘brmula (28) is an alternative hut equivalent 
representation of the power coefficient originally 
derived by Rice n . 

Results 

We illustrate the usefulness of TBIEM3D in 
computing the generalized radiation impedances and 
power transmission coefficients by considering a wide 
range of operating parameters and mode numbers. 

In all results the duct radius and length are one 
meter. Incident sound is generated by a circular array 
of either spinning or non-spinning axial dipole point 
sources located just aft of the duct trailing edge. 
Variations in impedance and power transmission due 
to duct length and source type were not investigated. 

Results are presented for nine cases. Three 
circumferential mode numbers, m = 0, 4, and 10, are 
considered. For each circumferential mode, we 
calculate results for three inflow Mach numbers, M = 
0.0, 0.2, and 0.6. All calculations were performed on a 
personal computer with a Pentium 133 processor with 
32 megabytes of RAM.. 

In figures 3-5, the real and imaginary parts of the 
generalized radiation impedances Coo * * and C°n 

are plotted as functions of k for the three inflow Mach 
numbers. Figures 6-8 contain graphs of Ch > Cn * and 
C 4 22 . The functions C\°n Cn* and Cn are displayed 
in figures 9-11. Vertical lines in these figures mark 
the infinite duct cut-off frequencies. 

The power transmission coefficients n£, U° n and 

n" are plotted as functions of k in figure 12 for each 
Mach number case. Figures 13 and 14 display similar 
results for the m = 4 and m = 10 cases, respectively. 

Discussion 

Several researchers have studied modal 

impedances and power transmission at cylindrical duct 
terminations (e.g., refs. 3-8). In references 3-5, 
impedance calculations were based on Rayleigh’s 
formula. Later research such as that in references 6-8 
employed more precise boundary element or Wiener- 
Hopf noise prediction methods. The highest 
circumferential mode order for which calculated results 
were presented was m = 3. 

The results shown here for m = 0 and m = 4 agree 
qualitatively with those earlier results. Specifically, 
the direct impedances, all approach 1 + Oi for 

increasing k . Most discontinuities in the impedance 
curves occur at or near the cut-on frequencies of higher 
order radial modes. Except for the M = 0.6 case, the 
associated power transmission coefficients approach 


unity fairly quickly once the mode is cut on. 

For M = 0.0 and 0.2 , the cross-modal impedances 
Cl, and C7: in figures 3, 4, 6, 7, 9, and 10, are 
nonzero for smaller values of k and slowly tend to zero 
for large k . For M = 0.6 , intermodal coupling is quite 
weak for all values of k as evidenced by the results in 
figures 5, 8, and 11, These observations reinforce the 
intuitive concept that cross-coupling between radial 
modes should be relatively weak except when one of 
the modes is just above its cut-on frequency. 

Trends in the computed results are most clearly 
identified by examining the power transmission 
coefficients in figures 12-14. For example, there is 
little change in the transmission coefficients for the m 
- 0 and 4 cases as M increases from 0.0 to 0.2. 
Figures 12c, 13c, and 14c indicate reduced modal 
transmission as higher order radial modes are cut on. 
This effect is enhanced for larger order circumferential 
modes. 

As the circumferential mode order m and inflow 
Mach number M increase, the generalized radiation 
impedance and power transmission curves become 
more oscillatory. The reasons for these phenomena are 
not fully understood and require further study. 

Johnston and Ogimoto (ref. 7), noted considerable 
influence of duct length on the radiation impedances. 
This is absent in our calculations because of the 
shortness of the duct examined. In reference 7, a duct 
length of 10 radii was used, while the current study 
considers a duct length of one radius. The short duct 
greatly weakens the organ pipe resonance phenomenon 
responsible for the duct length effects observed in 
reference 7. 

Concluding Remarks 

We have examined two aspects of duct 
propagation and radiation which can contribute to 
more efficient fan noise predictions. Rayleigh’s 
formula was expanded to include uniform axial inflow. 
Forward and aft radiation can be computed using 
equation (9). Comparisons with TBIEM3D results 
shows that Rayleigh’s formula with single radial mode 
input accurately predicts sound directivity and levels in 
the peak noise regions both forward and aft. 

The ducted fan noise prediction code TBIEM3D 
was used to calculate modal impedances and 
transmission coefficients for a wide range of operating 
parameters. Results were obtained for higher Mach 
numbers and circumferential mode orders than has 
been previously published. Due to their simplicity and 
ease of use, both of these tools should be of value to the 
engine nacelle designer. 
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Figure 1: Comparison Between Rayleigh s Figure 2: Comparison Between Rayleigh s 
Formula and TBIEM3D Formula and TBIEM3D 

M = 0.0 m =20 k = 25.0 M = 0.4 m = 20 k = 27.3 
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Figure 7: Generalized Radiation Impedance 

Figure 8: Generalized Radiation Impedance 

m = 4 M = 0.2 

m = 4 M = 0.6 

Resistance 

Resistance 

Reactance 

Reactance 
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Specific Acoustic Impedance (Re p 0 c) Specific Acoustic Impedance (Re p 0 c) Specific Acoustic Impedance (Re p 0 c) 



Figure 9: Generalized Radiation Impedance Figure 10: Generalized Radiation Impedance 
m = 10 M = 0.0 m = 10 M = 0.2 

Resistance Resistance 

Reactance Reactance 


12 

American Institute of Aeronautics and Astronautics 



Specific Acoustic Impedance (Re p 0 c) Specific Acoustic Impedance (Re p 0 c) Specific Acoustic Impedance (Re p 0 c) 
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a) M = 0.0 a) M = 0.0 





